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The Photoluminescent (PL) measurements on photons and the transport measurements on ex- 
, citons are the two types of independent and complementary detection tools to search for possible 

exciton superfluids in electron-hole semi-conductor bilayer systems. In fact, it was believed that the 
■ transport measurements can provide more direct evidences on superfluids than the spectroscopic 

_ measurements. It is important to establish the relations between the two kinds of measurements. 

In this paper, using quantum Heisenberg-Langevin equations, we establish such a connection by 
calculating various exciton correlation functions in the putative exciton superfluids. These cor- 
, _ relation functions include both normal and anomalous Greater, Lesser, Advanced, Retarded, and 

, Time-ordered exciton Green functions and also various two exciton correlation functions. We also 

evaluate the corresponding normal and anomalous spectral weights and the Keldysh distribution 
functions. We stress the violations of the fluctuation and dissipation theorem among these various 
exciton correlation functions in the non-equilibrium exciton superfluids. We also explore the input- 
output relations between various exciton correlation functions and those of emitted photons such as 
the angle resolved photon power spectrum, phase sensitive two mode squeezing spectrum and two 
photon correlations. Applications to possible superfluids in the exciton-polariton systems are also 
C/2 ■ mentioned. For a comparison, using conventional imaginary time formalism, we also calculate all the 

exciton correlation functions in an equilibrium dissipative exciton superfluid in the electron-electron 
coupled semi-conductor bilayers at the quantum Hall regime at the total flUing factor vt ~ 1- We 
stress the analogies and also important differences between the correlations functions in the two 
exciton superfluid systems. 
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INTRODUCTION 



An exciton is a bound state of an electron and a hole in the band structure of a solid. When excitons are sufSciently 
apart from each other, they may behave as bosons. Although Bose-Einstein condensation (BEC) of excitons was 
proposed more than 3 decades ago^, no exciton superfluid phase has been observed in any bulk solids yet. Recently, 
,—1 . degenerate exciton systems in quasi- two-dimensional semiconductor GaAs/AlGaAs electron- hole coupled bilayers 
\^ ' (EHBL)i have been produced by many experimental groups with photo-generated method** or gate- voltage generated 
\ method^i^ and electron-electron coupled bilayers at quantum Hall regime at the total filling factor vt = ^ (BLQH)^"— . 
■ It was widely believed that EHBL and BLQH are two of the most promising systems to observe BEC of excitons 
T-H ] among any solid state systems. Indeed, there have been extensive experimental searches for exciton superfluids in 
T-H . both systems^i^ii^. Now it was more or less established that an exciton superfluid subject to substantial dissipations 
^""^ ■ has been observed in BLQH systerrbiSrJ^. But there are still no convincing experimental evidences that the exciton 
^ . superfluids (ESF) have been formed in EHBL at the present experimental conditions^"—. In the photo-generated 
electron-hole bilayer (EHBL) samples^, various kinds of photoluminescence (PL) measurements were made by different 
groups^. More recently, taking advantage of the long lifetimes of the in-direct excitons, researchers started to be 
able to manipulate the excitons movements and perform various transport properties of the excitonsi^"— . In the 
undoped electron-hole bilayer (EHBL) samples^'S which is a heterostructures insulated-gate fleld effect transistors, 
separate gates can be connected to electron layer and hole layer, so the densities of electron and holes can be tuned 
independently by varying the gate voltages. Transport properties such as the Coulomb drag or counterflow can be 
performed in this experimental set-up. The PL signals, although weaker than those in the photo-generated samples, 
is still measurable. The transport measurements in^i^ii^""— and the PL measurements in^ are complementary to each 
other. It is important to search for exciton superfluid from both experimental methods. Any signatures of the exciton 
condensations should show up in both type of experiments. In fact, it was believed that the transport measurements 
can provide more direct evidences on superfluids than the spectroscopic measurements. Indeed, it is all the peculiar 
transport properties which proved the existence of superfluid in the liquid ^He at very low temperatures. 

There are previous theoretical works studying various photon emission spectra from the putative ESF—""—. I n^^d'^ , 
assuming the ESF has been formed in the EHBL, the authors studied the angle resolved photon spectrum, momentum 
distribution curve, energy distribution curve, two mode phase sensitive squeezing spectrum and two photon correlation 
functions from such an ESF. They found these photoluminescence (PL) display many unique and unusual features not 
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FIG. 1: The irreversible photon emission process in an non-equilibrium exciton superfluid. The Time reversal symmetry is 
broken. The photons are emitted and never flow back, so the arrows only go to one direction. The outgoing photons correlation 
functions are related to those of excitons by the input-output relations discussed in Sec. VI. As shown in Sec. IV and V, the 
fluctuation dissipation theorem (FDT) does not apply to this kind of quantum open systems. 
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FIG. 2: The reversible exchange processes between exciton superfluid and photon bath in an equilibrium exciton superfluid. 
The Time reversal symmetry is not broken. The photons are flowing back and forth, so the arrows go to both directions. As 
shown in Sec. IX, the fluctuation dissipation theorem (FDT) holds in this kind of equilibrium systems. The dissipative exciton 
superfluids in BLQH^'^"^^ is in this class with the photon bath replaced by gapless electron- hole excitations ( see-^- ). 



shared by any other atoniie or condensed matter systems. Observing all these salient features by possible future angle 
resolved power spectrum, phase sensitive homodyne experiment and HanburyBrown-Twiss type of experiment could 
lead to conclusive evidences of exciton superfluid in these systems. However, all the previous theoretical works focused 
on the the PL of the emitted photons, but various properties of the excitons themselves have not been addressed. 
There are also some previous theoretical worki^i^ on the Coulomb drag in electron-hole bilayer in the BCS side by 
using weak coupling mean field BCS theory. However, as analyzed irti^iii, for the experimental relevant density regime 
n ^ lO^'^ cm^^, the excitons are tightly bound pairs in real space, so all the experiments are in the BEC regime, so it 
remains interesting to study the exciton correlations in the strong coupling BEC limit. 

As mentioned above, the EHBL and BLQH are two of the most promising systems to observe BEC of excitons 
among any solid state systems. However, there are some crucial differences between the exciton superfluids in the two 
systems. The exciton -|- photon system in the EHBL is a typical quantum open system subject to a dissipative bath 
shown in the Fig.l. The photon emission process is a typical irreversible outgoing process breaking the time- reversal 
symmetry. However, the exciton superfluid in the BLQH at the total filling factor i/^ = 1 is a typical equilibrium 
system, also subject to a dissipative bath shown in the Fig. 2. There are reversible exchange process between the 
excitons and the dissipative gapless excitations. These processes preserve the time reversal symmetry. The relations 
and differences between the two kinds of systems have never been addressed before. As mentioned above also, 
the PL measurements and the transport measurements are two independent and complementary measurements to 
search for exciton superfiuids in the EHBL. The PL detects the properties of the emitted photons, while the transport 
measurements detect those of the excitons. The excitons and photons are coupled to each other through the interaction 
term in the Hamiltonian EqnlTJ So it is important to study the relations between the photon correlation functions 
and those of excitons which will establish relations between the two kinds of measurements. These relations are 
also crucial in bridging the interconnections between the optical communications and the electronic signal processing 
required at the integrated circuit devices. Here, we will address these outstanding open problems. 

Specifically, using quantum Heisenberg-Langevin equations, we will study various exciton correlation functions 
in the steady photon emission state. We will also investigate various input-output relations between these exciton 
correlation functions and those of emitting photons. The results achieved on exciton correlation functions are directly 
relevant to the transport properties of the excitons. They are also needed to evaluate the superfluid density and 
the critical velocity of a non-equilibrium exciton superfiuid in Fig.l and equilibrium exciton superfiuid in Fig. 2. The 
input-output relations can also be used to establish the connections between the transport measurements in^i^ii^rJ^ 
and the PL measurements We will also explore the analogies and differences between the exciton correlation 
functions in the non-equilibrium exciton superfluid systems shown in Fig.l and those in the equilibrium dissipative 
exciton superfluid systems shown in Fig. 2. Non-equilibrium physics remain poorly understood. So the results achieved 
in this manuscript in the context of non-equilibrium exciton superfluid may lead to some general concepts applicable 
to various other non-equilibrium problems, they could also shed some lights on the equilibrium dissipative exciton 
superfluid observed in BLQH systems--. 

In parallel to the experimental search for the exciton superfiuid in the EHBL, there are also extensive experimen- 
tal activities to search for exciton-polariton superfluid inside a planar micro-cavity. Although exciton condensation 
in a single quantum well ( SQW ) has not been observed so far, there are some evidences for the observation of 
Exciton-polariton condensation in SWQ enclosed inside a planar microcavity^"^. These evidences include macro- 
scopic occupation of the ground state, spectral and spatial narrowing, a peak at zero momentum in the momentum 
distribution and spontaneous linear polarization of the light emission and so on. One photon and two photon cor- 
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relation functions have been measured in the experiments respectively. We expect that the methods developed and 
results achieved in this paper should also apply to the exciton-polariton condensation inside a microcavity. 

The rest of the paper is organized as follows. In Sec. II, to be self-contained, we review the effective exciton-photon 
interaction Hamiltonian, the input state, the input-output relation and the exciton decay rate already studied ini^iii, 
but we will use a slightly different notations than^^'^^ which are consistent with the condensed matter notations 
in2^. Most importantly, we also write out explicitly the exciton operator which will be the starting point to calculate 
various exciton correlation functions in the following sections. In Sec. Ill, we calculate the one photon correlation 
functions, especially the one photon anomalous Green function which was not evaluated ini^iii. The central results 
of this paper are presented in Sec. IV-IX. We compute various normal one photon correlation functions such as the 
greater G^, lesser , retarded G^, advanced G^, the time-ordered and the Keldysh component G^ in Sec. IV, 
various anomalous one photon correlation functions such as the greater G^ , lesser , retarded G^, advanced G^, 
time-ordered G^ and the Keldysh component G^in Sec.V. Then we explore various Input-output relations between 
the photon correlation functions and those of excitons. In Sec. VII, we explore the non-trivial relations between the 
angle resolved photon power spectrum (ARPS) and the normal spectral weight and the normal distribution function, 
between the two mode squeezing spectrum and the anomalous spectral weight and anomalous distribution function 
at the experimentally relevant values. In Sec. VIII, we calculate the two exciton correlation functions and study their 
connections and differences than those of photons which can be measured by HanburyBrown-Twiss type experiments. 
These two-exciton correlations functions are important to transport measurements on excitons. In Sec. IX, using 
the imaginary time formalism, we study all the exciton correlations functions of an equilibrium dissipative exciton 
superfluid subject to the photon dissipation bath in Fig. 2. Then we compare with the results achieved in the previous 
sections on non-equilibrium exciton superfluid in Fig.l. We reach conclusions and possible future perspectives in the 
final section X. Throughout the paper, we did not consider the possible important effects of spins of electrons and 
holes which lead to the formation of the bright excitons with J = ±1 and the dark excitons^ with J = ±3/2, the 
effects of the trap, finite thickness and disorders in the sample. These facts will be briefly discussed in the conclusion. 



II. THE PHOTON-EXCITON INTERACTION, INITIAL STATE AND INPUT-OUT FORMALISM 

The total Hamiltonian of the exciton + photon is the sum of exciton part with the dipole-dipole interactions, the 
photon part and the coupling between the two parts Ht = Hex + Hph + Hint wher o^^'^'' : 

k kpq 

Hph = y^^u/kajak 

k 

H^nt = Y.^^g{k)ake'^^'b\ + h.c.] (1) 

k 

where A is the area of the EHBL, the exciton energy E'^ — k'^ /2M + Eg — Eb, the photon frequency Uk = Vg\J k^ + k"^ 
where Vg — cj -Jt with c the light speed in the vacuum and e ^ 12 the dielectric constant of GaAs^ k — (fc, fc^) is the 
3 dimensional momentum, Vdifi) is the dipole-dipole interaction between the excitonsi^, Vd(|r^ ^ d) = e^<P' j \ and 
= 0) = 22[^ where d is the interlayer distance leads to a capacitive term for the density fluctuationi^. The 

^ —1/2 ^ 

^(fc) ^ th\ ■ Dk X Lz is the coupling between the exciton and the photons where ekx is the photon polarization, Dk 
is the transition dipole moment and iz ^ oo is the normalization length along the z direction^^. In a frame rotating 
with the frequency /i, the Hamiltonian becomes 

k kpq 
k 

HL - Y.^i9{k)akb\ + h.c.] (2) 

fe 

In the following except in Sec. IX, all the calculations will be done in this rotating frame. 
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In the dilute limit, the average distance between the excitons is large, so the Vd is relatively weak, therefore we 
may apply the standard Bogoliubov approximation. In the Exciton superfiuid (ESF) phase, one can decompose the 
exciton operator into the condensation part and the quantum fluctuation part above the condensation: 



h = ^NSf:, + bj: (3) 

The linear term in bk is eliminated by setting the chemical potential ^ = Eq^ + nVd{0) — {Eg — Eb) + nVd{0). In a 
stationary state, the fi is kept fixed at this value. Upto the quadratic terms, the exciton Hamiltonian Hex is simplified 
to: 



k 



[(e,~ + Vd{k)n)b'^-, + i^^bW_. + h.c.)] (4) 

where the density of the condensate n ~ N/S. It can be diagonalized by Bogoliubov transformation: Hgf = E{0) + 
where E{k) = \J tf:[tj: + 2nVd{k)] is the excitation spectrum, (3^ = uj^bj: + v-j^'^ - is the Bogoliubov 



quasi-particle annihilation operators with the two coherence factors = ± -I. As fc — > 0, E{k) = u 

k k '2,E\kj 



shown in Fig. 3 where the velocity of the quasi-particle is it = yJnVd{Q)/M . At T = 0, the number of excitons out 
of the condensate is: n\T = 0) = ^ X]a;(^^%) — I (tlry^ which is the quantum depletion of the condensate due to 
the dipole-dipole interaction. One can decompose the interaction Hamiltonian Hint in EqnlT] into the coupling to the 
condensate part -fff„j = J2k [w{kz)VNak^ + h.c] and to the quasi-particle part Hf^^ — Y^^[ig{k)aiJb\ + h.c]. In the 

following, we study fc = and k ^ Q respectively. 

The initial state of the exciton-|- photons system in Fig.l is taken to be: 

\in) = \BEC)\(S)ph (5) 

Assuming a non-equilibrium steady exciton superfiuid state has been reached, the output field a^"*(a;) in the Fig.l 
is related to the input field by^^i^i: 

af\Lo) = [-l + i7,-G„(fc,c.-f*^)]al"M 

-«7fcGa(fc,w + iy)a"l(-a;), (6) 

where the normal and anomalous Green functions are^i^: 

LO + €j^ + nVd{k) 



G„(fc, w) 



E^(k) 



GJfc,..) = (7) 
^ ' u^-E^{k) ^ ' 



where the cj = cj/j — /i. The exciton decay rate in the two Green functions are^^: 



|2 



7fc = 27r£ij:(^) = ^)| (8) 

which is independent of L^, so is an experimentally measurable quantity. Just from the rotational invariance, we can 
conclude that 7^ ~ const. + as fc — >■ as shown in Fig. 3. 
From Ref,-^, one can also find the exciton operator—: 

- Ga(fc,.. + z^)a3(-c.)] (9) 

Note that the input-output formalism Eqns i6l9] stands for a non-equilibrium dynamic process. It takes Eqnl5]as the 
initial state, then solving the coupled Heisenberg equations of the photon and exciton operators. The steady state 
of the excitons has been reached, the photons are scattering from such an exciton steady state, so the Eqn[S] can be 
considered as the scattering S matrix shown in the Fig.l. 

In the following sections, we will first use the Eqn|5] to calculate the one photon normal and anomalous Green 
functions, then we will use the EqnjS] to evaluate all the 5 normal and anomalous exciton correlation functions: the 
greater G^, lesser G^, retarded , advanced G^, the time-ordered G^ and the Keldysh component G^ . 
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III. ONE PHOTON NORMAL AND ANOMALOUS CORRELATION FUNCTIONS 

From EqnjSl one can compute the photon umber spectrunii^: 



(10) 

n2(a;)+7|£;2(fc) 



By using the Fourier transformation to S'i(fc, w), one can get G±{k, r) = G'i(fc, r) 



where n{uj) = u? - E'^{k) + 



which, obviously, only depends on the anomalous Green function Ga{k,uj). 

The normalized first order correlation function g±{k, r) = 5i(fc, r) = Gi(k, t)/Gi {k,T = 0). When t > 0, it is given 
by: 

g^{k,T) = e-^^[cos(i;(A^)T) + sm{E{k)T)] (12) 

2E{k) 

It turns out that the first order correlation function is independent of the relation between E{k) and 7jr/2. 
Similarly, one can find: 

(a°-"*(a.)af = | - 1 + 7,,G„(fc,a. + z^^)]' (13) 

which, in sharp contrast to EqnfTTl only depends on the normal Green function G„(fc, w). 

The power spectrum experiments^^ only measure the normal ordered one photon correlation function Eqn llll The 
anti- normal ordered one photon correlation function Eanll3lmavbe measured by a transmission spectrum with a weak 
external pumping. 

One can also calculate the anomalous one photon Green function: 

{aZ-\uj)a1-i{-u)) = [-l+7-G„(fc,.. + z^)] 

X j^Ga{k,~u + i^) (14) 

which depends on both the normal and the anomalous Green functions. 
By changing k — —k,u; — > —uj in Eqn ll4l one can get: 

(a°_|(-w)a°-"*(..)) = [-l+j^G,,{k,^u + t^)] 

X j^Ga{k,UJ + t^) (15) 
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From the explicit forms of Gn,Ga listed in Eqnl?! one can observe that the right hand side of EqnlT4]is a even 
function of w, so: 

(a2"*(c.)a°|(-c.)) = (a°|(-c.)a2"*(c.)) (16) 
This identity is expected because both the input and output fields obey the Bose commutation relations^: 

[a^{t),aZ^,Ht')] = [af (t),a°r^(t')] - hl^t - t') 

[af{t),a^{t')] = [ai:^\t),alfit')] = (17) 
so that (ae"*(t)a°"i(0))„ = (a°"i (0)a™*(i))jn whose Fourier transform leads to EqnUl 

k —k —k k 

One can also find the Fourier transforms of EqnlT4] and [151 

Fi{r) = I ge-(a°_|(-a.)ae"*(a.)) 

F2ir) - J ^e-^(a2"*(c.)a°_|(-c.)) (18) 

Using EqnfTHl one can see that Fi{t) = F2{t). Then the normalized anomalous one photon correlation function in 
the time domain is given by: 

Mr) = fi(r)M(0)^- '''(^^ + ^ e-^- 

nVd(k) 

~iE{k)T piE{k)T 

X H lr. ... +4 J.. (19) 

2 



which was used in calculating the two photon correlation functions in^^. 

Measuring these photon anomalous Green functions are through the phase sensitive homodyne measurements dis- 
cussed in thei^. 



IV. ONE EXCITON NORMAL CORRELATION FUNCTIONS 

From EqnlSl one can compute the normal exciton correlation function: 

iG<{k,u) = {bl{u)bj:{u)),r,^j~^\Ga{k,co + z^)\^ 



^^^^ 



n^{iu)+-fiE'^{k) 



where ft{uj) — uP' ~ E'^{k) + 72/4. Note that it only depends on the anomalous Green function. Obviously, the 



normalized first order exciton correlation function 



guik, t) = Gib(fc, r)/Gib(fc, r = 0) = gi{k, r) (21) 



is also given by EgnlT^ and shown in Fig. 4. 

From Eqnl9l one can compute the exciton correlation function: 



zG>(fc,w) = {b-,{u)bl{u)),„^^f:\G,,{k,u + t^)\^ 

_ 7g[(^ + 6g + ny,(fc))^ + (7g/2)^] ^22) 



n^{Lu)+-fiE^{k) 

which only depends on the normal Green function, in sharp contrast to Ean l20l 
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r(ns) 

FIG. 4: The one exciton correlation function of Ean l21l at the experimental relevant valued nVd{k) = 5fj,eV, 7j;/2 = IfieV at 
the quasi-particle energies E{k) — 2fj,eV. The same values of the nVd(k),^%/2 will also be used^^ in Fig. 5-10. 



From Eqnl^and Eqiil^Ol one can find the normal spectral weight 

Pn{k,uj) ^ {b^{uj)bi{uj))in - {bi{uj)bj:{uj))in 

= 7fe[|G„(fc,c. + z^)|2-|G,(fc,w + *^)n (23) 
Using the Eqnl?! one can manipulates the last equation in [23] into a surprisingly simple result: 

Pn{k,uj) = i[Gnik,uj + iy) - G„(fc,w - iy)] (24) 
One can also calculate the Retarded and Advanced normal Green function: 

iG^{k,t) = e{t)[b-,{t),blm 

iGt{k,t) = -e{-t)[b-^{t),b\m (25) 
Their Fourier transforms lead to their corresponding normal spectral weight representations: 

duj' pnik^uj') 



G^{k,oj) 



27r LJ — uj' + irj 



G-(fc>)=/|^^4^ (26) 
J Ztt u! — ui' — trj 

Plugging the Eqnl21]into the above equations leads to: 

G^{k,Uj)^Gn{k,UJ + t^) 

G^ik,u;) = G,,{k,u:-t^) (27) 

whose physical pictures are clear: the Retarded or Advanced normal Green functions can be achieved simply by just 
putting the decay rate ±^ in the non-dissipative ones in Eqnl?! 
From Eqns l23l24l27l we can check the identity: 

C'n ^ = G^ ^ (28) 

as expected. 

We can also compute the Time ordered exciton Green function 

iGliKt) ^ {Tbj:{t)blm 
- 0{t){b~^{t)b\m + e{-t){b\{Q)b~^{t)) (29) 
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Its Fourier transform leads to: 

j 2^^^73t;^"^7^i:;t7^J ^^^^ 

Substituting Egnl^ into the above equation and paying special attentions to the pole structures in the EgnlSOl 
one can show that 



1{k,uj) = [l-i7j:G„(fc,w + iy)]G„(fc,cj-Zy) 

= [l-i^^^G^{k,uj)]Gt{k,Lo) (31) 



Finally, from Eqnl5^ and Eqnl^Dl we can easily get the Keldysh component of the Green function: 

iG^{k,u:)^i[Gl{k,u^) + G<{k,uj)\ 
^ liii^ + H + nVa{k)f + {nVa{k)f + (7g/2)^] 
02(a;) +7Hi;2(fc) 

which is shown in the Fig. 6. It is neither a even nor an odd function. 

It is well known that for an equilibrium system, the and G^ are related by the Fluctuation-Dissipation Theorem 
(FDT) instead of being independent. At T = 0, the FDT takes the form 

i[Gl{k,u) + G<{k,uj)]^isgnu:[Gl{k,u)^G<{k,u)] (33) 

However, for a non-equilibrium system with the initial input state EqnlSl the G^ in Eqn l22l and the G^ Eanl20l 
violate the FDT EqnlSl EqnlMl and |30] lead to G^, G^ and G^ respectively. This is an salient feature expected 
for any non-equilibrium system. The initial input state EqnIS]is not even an eigen-state, let alone the ground state of 
the total Hamiltonian EqnlT] It is this fact which leads to the dynamic photon scattering process encoded in Eqnl6] 
andO The spectral weight p„(/e,a;) in EqnUSland the Keldysh component G^(fc,w) in Eqnl32] automatically follow. 
Detailed comparisons with the corresponding equilibrium system will be given in Sect. IX and the conclusion section. 

V. ONE EXCITON ANOMALOUS CORRELATION FUNCTIONS 

From Eqnini one can also compute the abnormal exciton correlation function iG^{k,uj) — {h^{uj)h 



(34) 



iG>{k,uj) = 7jG„(fc, w + iy)Ga(/i:, -w + i—) 

_ + eg + nVd{k))nVd{k) 

Q'^{uj)+-fiE^{k) 

which depends on both the normal and the anomalous Green function. 

By changing k — > —k,uj ~uj in Eqn l341 one can get iG'^{k,Lu) = {h_^{—Lj)h^{uj))in: 

^G<{k,u;)^^J:Gn{k,^u; + ^^)Ga{k,iu + ^^) (35) 
From Eqn l34l and 1351 one find the abnormal spectral weight 

p„(fc,c.) = z[G>(fc,w)-G<(fc,L.)] (36) 
When using Eqn ll4[ 1151 and [TBI one can manipulate it into a surprisingly simple result: 

Pa{k,Oj)^l[G>{k,Lj)~G<ik,io)] 

^i[Ga{k,UJ + t^)-Ga{k,UJ-l^)] (37) 

In contrast to Eqn[T6]for photons, G-^{k,uj) 7^ G^{k,uj), so pa{k,uj) 7^ 0. This is expected, because {b^{t)b_^{0)) 7^ 
(^-£(0)%(i))- 
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One can also find the Fourier transforms of Eanl35]andl34l 



F2,{r) = j ^e-zG>(fc,a.) (38) 

In contrast to the corresponding quantities for photons Fi{t) = F2{t) hsted in EgnfTSl Fh,{t) ^ F2b{T). In fact, one 
can see that Fi,(t) = ^^i(t)/7,^, F2,(t) = Fu{t) - znFrf(fc)e-^- 

Then the normalized anomalous one exciton correlation function in the time domain is given by: 

fibir) = Fib(T)/Gib(0)-/i(T) 



f2b{r) = F2b{T)/Gu{0) = /i,(r) - t2 



nVdik) 



m 



which will be useful in calculating two exciton correlation functions in the Sect. VIII. 
One can also calculate the Retarded and Advanced abnormal Green function: 

Their Fourier transforms lead to their corresponding spectral weight representations: 



iG^iKt) = -9{~t)[b~^{t),b_~^{0)] (40) 



Ga {k,^) = / -X 7— — 

J Ztt Lo — ui + trj 

Gtik,.)=[^J^^^ (41) 

J ZTT U! — LO' — IT] 

Plugging the Eqn|37]into the above equation leads to: 

Gi{k,uj) = Ga{k,uo-i^) (42) 

whose physical pictures are clear: the Retarded or Advanced anomalous Green functions can be achieved simply by 
just putting the decay rate ±^ in the non-dissipative ones in EqnjTl 
From EqusinZl we can check the identity: 

Ga ~ Ga — G^ — G^ (43) 

as expected. 

We can compute the Time ordered anomalous exciton Green function: 

zGlik,t) = {Tbj:{t)b_^,iO)) 
= 0mbj:{t)b_j:{O)) + ei-t){b_^^{0)bj:{t)) (44) 

Its Fourier transform lead to: 

r,r . [dJ^ G>(k,J) G<{k,u') ^ 

Ga{k,i^j)^ —[— : — J (45) 

J Am ui — ui — ir\ u! — uj + ti] 

Substituting Ean l34l55] into the above equation and paying special attentions to the pole structures in the Eqn l45[ 
one can show that 

Gl{k,U;) - [l^ij^^Gnik,LU + l^)]Ga{k,U-t^) 

= G^ik,uj) + G>ik,io) 

= G«(fc,^) + G<(fc,a;) (46) 
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Finally, from Eqnl22] and Eqnl20l we can easily get the Keldysh component of the anomalous Green function: 

iG^{k,io)=i[G>{k,uj) + G<ik,io)] 
27g(^^ + eg + nVd{k))nVa{k) 
n'^{u)+jiE^{k) 



(47) 



which is a complex even function. Its real and imaginary part are shown in the Fig. 8 and Fig. 9 respectively. 

Similar to the normal exciton correlation functions discussed in the previous section, as expected for an non- 
equilibrium system, the G^ in EgnlMl and G^ Eqnl25]are two independent Green functions, so violates the FDT for 
anomalous Green function: 

i[G>{k,u;) + G<{k,Lo)] = isgnu:[G>{k,u:) - G<(fc,c.)] (48) 

Eqn|3T]and|3S]lead to G^, G^ and G^ respectively. The spectral weight /9,i(fc, lo) in EqnI27]and the Keldysh component 
G^{k,uj) in Eanl47l automaticallv follow. Detailed comparisons with the corresponding equilibrium system will be 
given in Sect. IX and the conclusion section. 



VI. INPUT-OUTPUT RELATIONS BETWEEN PHOTONS AND EXCITONS 

From the exciton condensation Eqnl3]at fc = and the emitted photon number at fc = Eqn.4 in Ref llOl ( which is 
the same as the photon condensation Eanl62l ]. one can immediately see the relation between the two: 

Nph = Nbjo (49) 

where the 70 is the exciton decay rate at fc = 0. This relation relates the condensation of the emitted photon to that 
of excitons. 

From Eqn llOI201 one can see that 

(af ^(^)a2"*(^))„ = 7fc(&t(^)&£(^))» (50) 

which relates the emitted photon spectrum to the internal normal ordered exciton correlations. A similar relation 
was discussed in the context of optical cavities^^iSSi. Here it is in the spontaneous photon emission from an exciton 
superfluid in the absence of any cavities. This relation shows that the angle resolved photon power spectrum studied 
ini^ can reflect precisely the normal ordered normal exciton correlation function iG^ . 

However, one can also see that due to the extra first term 1 in Eqn ll31 there is no such simple input-out relation 
between (ae"*(a;)aT*^(w)) and {bj:{uj)bt{uj)) listed in EqnUl 

Now we will try to explore the input-output relation between the anomalous photon correlation function in Eanll4l 
and that of the excitons in Eanl34l Obviously, they are not simply related as the Egn lSOl However, by comparing 
Ean ll4l with Ean l461 one can establish the input-output relation between the anomalous photon correlation function 
with the time-ordered correlation function of the excitons: 

(a°-«*(a.)a°|(-a.)) = J^,^Gl{k,u;) (51) 

which relates the emitted photon two mode squeezing spectrum to the internal time-ordered exciton correlations. A 
similar relation was discussed in the context of optical cavities^ Here it is in the spontaneous photon emission 
from an exciton superfluid in the absence of any cavities. This relation shows that the phase sensitive homedyne 
measurements to measure the two mode squeezing spectruroii can reflect precisely the Time ordered anomalous 
exciton correlation function iG'^. It is instructive to compare EgnlSlI which involves the Time-ordered with EgnlSOl 
which involves just normal ordered. 



VII. THE QUASI-PARTICLE SPECTRA AND DISTRIBUTION FUNCTIONS IN A 
NON-EQUILIBRIUM STATIONARY EXCITON SUPERFLUID 

There are three kinds of Photoluminescence measurements. One kind is the angle resolved power spectrum (ARPS) 
shown in Fig. 3 ini^. The second kind is the two mode squeezing spectrum shown in Fig.lb-3a irtii. In any non- 
equilibrium system, one need both the spectral weight Pn/aik, to) — i[G,^/ (fc, w) — G^,(k, w)] and the Keldysh Green 
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FIG. 5: The normal spectral weight Ean l52l at different quasi-particle energies at E{k) < 7g/2: E{k) — 0.05, 0.5 and at 
E{k) > 7g/2: E{k) = 2,5,8^61/. 



function component G^^^{k,Lu) = G^^^{k,uj) + G'^^^{k,uj) to characterize the spectral weight and the distribution 
function respectively. In this section, we will explore the relations between the two kinds of PL with the exciton spectral 
weights and distribution functions. Due to the lack of FDT in the non-equilibrium exciton superfluid system, these 
relations are non-trivial, so need to be studied in details. The third kind is the two photon correlations measurement 
shown in Fig. 3b in-i^. In the next section, we will explore its relation with the two exciton correlation functions. 



A. Connections among the ARPS, the exciton normal spectral weight and distribution function 
From EqnlMl we can get the normal spectral weight pn{k,uj) = —2ImG!^{k,uj): 



'■{u:-E{k)r + {jf:/2f 

2 7£ 

{u^ + E{k)Y + {^j:/2Y 



-I , . .... (52) 



where ""I ~ — 1- Obviously, pn{k,uj) satisfies the sum rule: 

^p„(fc, uj) ^ul-vl = l (53) 

as expected from its definition. 

It is shown in Fig. 5 for several different quasi-particle energies E{k). At E{k) — 0, the p„(fc = 0, w) = Lj'i+'(^^j2y^ ^ ^ 
is a single Lorentizan centered at w = with width 7j:/2. At E{k) > 0, it consists of two Lorentizan centered at 
zLE{k) with width 7j:/2 with the spectral weight u| and —vi respectively. Only when E{k) ^ 7^/2 in Fig. 5, namely, 
when fc >• fc* in the Fig. 3, the two Lorentizans are well separated. Note that oj — 0) = ~ — — — > is large 

when E{k) < 7j:/2. The positive value at w = can be seen in Fig. 5 when E(k) < jj:/2. The normal distribution 
function iG^ is shown in Fig. 6. 

What the first class of PL experiment measured is the angle resolved power spectrum (ARPS) shown in Fig. 3 in — . 
The excitation spectrum and the distribution function of the exciton superfluid system to be probed by transport 
measurements is given by Eqnl52] and [32] respectively. So it is important to explore the relation between what 
is measured by the PL and the system's properties to be probed by transport measurements. Such a relation is 
given by the input-output relation Egn lSOl Ean l20l and Eanl52l One can see qualitatively the same physics with 
slightly quantitatively different numbers. When k < k*, the Fig. 5 shows that the quasi-particle is not well defined. 
Correspondingly, the ARPS in the Fig. 3 in-^*' has just one peak pinned at uj^ — p with the width 7fc. The MDC has 
large values at fc < fc*. The EDC has large values at w < 7^. From the Fig. 6, one can see that the distribution function 
has just one peak pinned at cok — p with the width 7^. When k ^ k* , the Fig. 5 shows that the two Lorentizans 
centered at ±E{k) with the width 7^:/2 are well separated with two different spectral weights m| and — w| respectively. 

The quasi-particles at zLE{k) are well defined. Correspondingly, the ARPS in the Fig.l iniS has two well defined 
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FIG. 6: The normal Keldysh component iG^ at difTerent quasi-particle energies E{k) = 0.05, 0.5, 2, 5, 8/ieV. 



symmetric quasi-particles peaks at 0;^ = /i ± [i?^(fc) — 7|/4]"'^/^ with the width 7fe/2. From the Fig. 6, one can see that 

the distribution function spHt into two asymmetric peaks at uji^ — fi ± [E^{k) — ji/i]^/^ with the width 7fc/2. It is 

instructive to observe some differences between the Pn{k, oj) shown in the Fig. 5 and the ARPS shown in the Fig. 3 in^^..: 
(1) There are two symmetry peaks in the latter, but the two peaks in the former are as-symmetric with two different 
spectral weights ui and — subject to the constraint u| — = 1. (2) There are very slight differences at the two 

peak positions, the latter are at ±[E^{k) — 7,-/4]^^^, while, the former are at ±E{k). The analogies and differences 



between the iG^ shown in the Fig. 6 and the ARPS shown in the Fig. 3 ini^ can be similarly discussed. 



B. Connections among the two mode squeezing spectrum, the exciton abnormal spectral weight and 

distribution function 



Similarly, from Eqnl371 we can get the abnormal spectral weight pa{k,uj) — —2ImG^{k,uj): 

Ik 



iLo~E{k)r + ijj:/2)^ 
7fc 1 



{u; + E{k)r + ijf:/2r 



(54) 



where Uj:v^ 



nVdjk) 
2E(k) ' 



It is an odd function of ui and satisfies the sum rule: 



/ 



du) 

27T 



paik,Uj) = 



(55) 



as expected from its definition. 

It is shown in Fig. 7 for several different quasi-particle energies E{k). It also consists of two Lorentizan centered 
at ±E{k) with width 7j:/2 with the spectral weights ^uj:vj: respectively. Only when E{k) 3> 7j:/2, namely, when 

k ^ k* in the Fig. 3, the two Lorentizans are well separated. The pa{k,uj) shown in Fig. 7 is an odd function of lj. 
From Ean l47[ one can see that the anomalous distribution function iG^ is a complex even function, its real and 
imaginary parts are shown in Fig. 8 and Fig. 9 respectively. Because the ratio of the real part over the imaginary part 
is independent of w, so they have similar shape. 

What the second class of PL experiment measured is two mode squeezing spectrum shown in Fig.lb-3a in^^. The 
excitation spectrum and distribution function of the exciton superfluid system to be probed by transport measurements 
are given by Eqnl54]and|47l Such a relation is given by the input-output relation EqulSU Eqn|35]and Eqni54l Eqnl47l 
One can also see qualitatively the same physics with slightly quantitatively different numbers. When k < k* , the Fig. 7 
shows that the quasi-particle is not well defined. Correspondingly, the two mode squeezing spectrum in the Fig. lb inii 
has just one peak centered around ujk = p, with the width 5i(k) which depends on the dipole-dipole interaction nVd{k). 
The corresponding squeezing angle is always positive in all the frequency range as shown in Fig. 4a inii. From the Fig. 8 
and 9, one can see that the distribution function has just one peak pinned at ujk — p with the width 7fc. When k ^ k* , 
the Fig. 7 shows that the two Lorentizans centered at ±E{k) with the width are well separated with the spectral 
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FIG. 7: The anomalous spectral weight Ean l54l at different quasi-particle energies E{k) = 0.05, 0.5, 2, 5, 8fieV. 
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FIG. 8: The real part of the anomalous Keldysh component iG^ at different quasi-particle energies E{k) — 0.05, 0.5, 2, 5, SfieV. 

weights =pugug respectively. The quasi-particles at ±E{k) are well defined. Correspondingly, the two mode squeezing 
spectrum in the Fig. 3 in^- has two well defined quasi-particles peaks at = /i ± [E^{k) — ^2/4]^/^ with the width 
(52 (fc) which depends on the dipole-dipole interaction. The corresponding squeezing angle becomes negative at small 
\ujk — At| and increase to be positive at large \u}k — vanishes at the resonant condition \ujk — m| = — 72/4]^/^ 

as shown in Fig. 4a in~^. From the Fig. 8 and 9, one can see that the distribution function split into two symmetric 
peaks at Wfc = ± [E^{k) — ^"i/A^^I"^ with the width 7/c/2. It is instructive to observe some differences between the 

Fig. 7 and the Fig.lb-3a in~i: (1) The widths in the former also depends on the the dipole-dipole interaction nVd{k). 
While the width in the former is just 7fe. (2) There are very slight differences at the two peak positions, the latter are 
at ±[£'2(fc) - 72/4] where the squeezing angle also vanishes, while, the former are at zLE{k). The analogies and 

differences between the iG^ shown in the Fig. 8 and 9 and the Fig.lb-3a io^ can be similarly discussed. 

Recently, the elementary excitation spectrum of exciton-polariton inside a micro-cavity was measured^, and was 
found to be very similar to that in a '^He superfluid except in a small regime near k = 0. We believe this observation 
on the anomaly near fc = is precisely due to the excitation spectrum in a non-equilibrium stationary superfluid 
shown in Fig. 5 and Fig. 7. 

VIII. TWO EXCITON CORRELATION FUNCTIONS 

In this section, we will explore the relations between the two photon correlations measurements and the two exciton 
correlation functions. Similar to two photon correlation functions studied ir>ii, one can also define the two exciton 
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correlations functions: 



92b iV 



92b iV = 



bUt)bUt + T)b-,{t + r)bf:{t) 



|Gi.(0)|^ 
r)bj:it + T)b^ Jt)b_dt) 



|Gi.(0)r 



(56) 



where the Gibir) = (^^(i + ''')bj:{t))in is the single exciton correlation function in Egnl^ and shown in the Fig. 4. In 
order to directly contrast with the two photon correlation functions calculated in^- , we first consider the orderings of 
the exciton opertaors in EqnlSB] The other orders will be discussed near to the end of this section. 

Due to the input-output relation Egn lSOl one can see immediately that (?2?(''") ~ 9^2 i''') 'v^^hich is given in Eqn.l2 
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However, due to the highly non-trivial input-output relation Egn lSll the g^^\T) could be even qualitatively 



differ from g2 (t) given in Eqn.l2 in^^. This is indeed the case as shown below. 
By using EqnOand the commutation relations in the frequency space [a'-" 



one can find 



0) + Fu{t) X F*,{t) 



(57) 



where the -Fifc(r) and F2b{T) are given in Ean l39l In fact, one can physically understand Eqn l57l bv noting that the 
first term is due to the equal time single exciton correlation function Gib(T = 0) = {i>^_^^{t)b^^.{t)) in the Eqnl2n 
while the second term is due to the exciton anomalous Green function Eqn|31] which encodes the exciton phase 
correlations. Because Fib{T) ^ -Pjh (''■), so the G^2b {'^) must be complex. Then the normalized two exciton correlation 
9'^^''\t) = 1 + fib{T) X /2f,(T) where the fu{T) and /26(t) are given by EqnOH so 



(±fc)/ ^ 



92 



{r)-^2{ 



E^k) + 



nVd{k) 



^^, sm{E{k)T) 
E{k) 



(58) 



The HanburyBrown-Twiss type of experiments are designed to measure the two photon correlation function shown 
in Fig. 3b in^^. The two exciton correlation function of the exciton superfluid system is given by Eqn l58l Although 

the normalized two photon correlation function Q^^^t) ( shown in Eqn.l2 and Fig. 3b inil ) is always positive, due 

to the complex second term in EgniSSj the g^^\T) remains complex. We can separately evaluate both the real part 
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FIG. 10: Real and Imaginary part of the two exciton correlation functions Ean l59l for E{k) — 2^eV . 



and the imaginary part of the difference of the two: 



^eg^f (r) - 9i^'\r)=-2 



[sin^ E{k)T - 



inVd{k)y 
7fc 



{nVd{k)y 
[sm.2E(k)T + 



■iiin2E(k)T] 
iEik) 



7fc 
Eik) 



sin^ E{k)T] 



(59) 



which is drawn in the Fig. 10. 

When comparing the two exciton correlation function with the single exciton correlation function in Eqn l21l and 
shown in the Fig. 4, one can see that the envelop function in the former ( latter ) decays with the 7^: ( 7j;/2 ), the 

oscillation frequency in the former ( latter ) is 2E{k) ( E{k) ). 

One can see why g''2b {''') in Eqnl56l must be positive definite. Define A'^ = 6t(i)6t(t + r),A = b^{t + T)b^{t), 
so the numerator can be written as H = A'' A which, of course, is a positive definite Hermitian operator. Similar 
arguments can be used used to show the two photon correlation functions defined in^^ must be positive definite. In 
fact, one can calculate various two exciton correlation functions by shifting the orders of the exciton operators or 



put Time-ordered inside the EqnlSSl For example, one can consider (^bt{t + t)6^ j.{t)b_j:{t)b^{t 



if one define 



A^ = bi{t + T)b^_^{t),A ^ b_j:{t)bj:Xt + r), the argument can be written as H — A'^A, so it must be positive definite. 



Similar procedures as above lead to its value 1 + |/i6(t)|'^ — g2^^\T). One can show that the other two orderings lead 

to 1 + |/26(t)P which is positive definite and 1 + fli^ir) x /26(t) — 9*2\^'^\t) which is complex conjugate to Eqnl58l 
and shown in the Fig. 10 by adding a minus sign to the imaginary part. 

Various exciton correlation functions discussed in this section are needed to evaluate the transport properties of 
excitons such the Coulomb drag and the counterflow resistances. 



IX. EXCITON CORRELATION FUNCTIONS IN AN EQUILIBRIUM DISSIPATIVE EXCITON 

SUPERFLUID 



In this section, we will study the Exciton correlation functions in an dissipative equilibrium exciton superfluid in 
Fig. 2 and then compare with those in the stationary non-equilibrium exciton superfluid in Fig.l studied in all the 
previous sections. 

As shown in the Sec. II, in the rotating frame with the frequency /i, the quantum open system in Fig.l has the 
Hamiltonian EqnI5]and the effective quadratic Hamiltonian Eqn|31 In the lab frame, the equilibrium system in Fig. 2 
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have the same Hamiltonian. Staring from Eqnl2]and effective quadratic Hamiltonian Eqn|4l we will discuss k = and 
k respectively. 

At k — 0, the Hamiltonian for the photons are: 

Ho = y^^{vg\kz\ - ti)a\k^)a{k^) 

i^VN{g{k,)a~ g*{k,)a^) (60) 



Obviously, it can be rewritten as 

Ho = Y^{vg%\ - pi){a\k,) ~ {a{k,))*){a{k,) ~ {a{k,))) + const. 



I (I. w iajk.WN 

{a{kz)) = — rrn (61) 

where {a{kz)) is the expectation value of the photon annihilation operator in the ground state. Because the coupling 

constant gikz) ^ Lz — > in the — oo limit, so a proper regulization is needed to extract the physics embedded 
in this expectation value. This was done ini^. Following the regulization procedure ini^, we can find the average 
photon number in the ground state 

n^,^=N-fo5{vg\k,^-fi) (62) 

where the 70 = \g\'^D is the exciton decay rate at fc = where the D = Lz/vg is the photon density of states at 
k = 0. One can see that the 70 is independent of Lz and is an experimentally measurable quantitj»i^. The photon 
condensation EquIHUis the same as the emitted photon number at A: = Eqn.4 in ReffTUl 

Now we will calculate the exciton correlations at fc 7^ in the ground state of the interacting system. Obviously, 
the exact ground state of the interacting systems |G) in Fig. 2 is not the same as the initial input state in Eqnl5l 

\G) ^ \in) = \BEC)\Q)ph (63) 

Even we do not know the explicit form of \G), we are still able to calculate all the ground state correlation functions 
by the imaginary time path-integral method which is particularly suitable to calculate correlation functions in an 
equilibrium systemSfi. 

In the imaginary time path integral representation of the Hamiltonian Eqn[T] 



^ = dT[bl{T)drbj:iT) + aliT)drakiT)] 



k 

+ / dT{Hsf +Hph + H,nt) (64) 
"'0 

Because the total action Eqn|64] is quadratic, integrating out the photons out exactly leads to the effective action 
for the excitons only: 



^ex = ^^[&t(ia;„)(-iw„ + E(fc,iw„))6j:(itj„) 



+ {^^^bl{^UJr,)bl.{^^u^) + h.c.)] (65) 



where the S(fc, iujn) is the self energy due to the integration of photons: 



S](fc,z^„)^^ . 'gf^l' (66) 

^ lUJn - [UJk - Ai) 
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where the sum is over continuous modes of photons at a given in-plane momentum k. 

One can easily find the imaginary time ordered T normal and anomalous exciton correlation functions^: 



(iw„-S(fc,iw„))2-S2(fc) 



Glik,-) = ^ (67) 

(*W„-S(fc,*C^„))2-i?2(fc) 

where the average is with respect to the exact ground state of the system. Obviously, the initial input state EqnlS] is 
not even an eigen-state of the system, let alone the ground state. This is the main difference between the equilibrium 
correlation functions calculated here and those calculated in the previous sections. It is the goal of this section to 
contrast some analogies and especially, the crucial differences between this two sets of correlation functions. 

After doing the analytic continuation ia;„ — > w + 177 in Eqn 1661 and replacing the sum by an integral over the density 
of state J2k, ~^ I d-^kD^{u!k), one can see: 

f Dr:(wfe)|g(fc)P 

S](fc,c^ + iv)^Pj dujk J_ ^'^ + (68) 

where the P stands for the principle part of the integral and the exciton decay rate 7^ is: 

7fc = 27r y dwfcL)^(wfc)|g(fc)|2(5(w - (wfc - //)) 

^ 2^Z?j(M)|5fcK=/i)|' (69) 

where we made the Markov approximation used iiti^iii. The real part causes a small energy shift which, in fact, 
vanishes identically within the Markov approximation ini^iii. Substituting Eqnsl68l and [69l into Eqn l67l leads to: 

g^{k,iu) - Gnik,^ + ^Y^^^n*ik,^) 

g^{k,uj) = Ga{k,u: + i^) = g^*{k,uj) (70) 

which are identical to Eqn l27l and Ean l42l At first sight, this result may look surprising. However, one can understand 
it better by comparing the equilibrium dissipative exciton superfluid in Fig. 2 discussed in this section with the non- 
equilibrium exciton superfluid system in Fig.l discussed in the previous sections. Both systems are described by the 
same Hamiltonian EqnlT] and 01 so should have the same excitation spectrum. This leads to EquITD] from a very 
general physical picture. However, what distinguishes the two systems is the initial conditions. Only the distribution 
function, namely, the Keldysh component depends on the initial conditions. While the retarded and advanced 
Green functions G^ and G"^ depends only on the Hamiltonian. The distribution function for a bosonic equilibrium 
system is simply given by the Bose distribution function, while that for an non-equilibrium system is determined by 
the initial conditions, namely the input state EqnjSl 

The greater and lesser normal Green functions at T = can also be found from the Fluctuation and Dissipation 
Theorem (FDT): 

ig>{k,u) = e{uj)pn{k,uj) 

ig<{k,uj) - -0{-uo)p^{k,uj) (71) 

Obviously, they differ from the corresponding Green functions in the non-equilibrium case Eqn l22l2ni Of course, 
we still have the identity: 

g^ - gt - gl - g^ (72) 

The real time-ordered T normal Green function can be found: 



gl{Kuj) = g^{k,u)e{u) + g:^{k,u;)9{~u) 

~ Gn{k,uj + i-^sgnoj) 

Obviously, it differs from the corresponding Green function in the non-equilibrium case Eanl30l 



Gn{k,uj + i-^sgnoj) (73) 
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From EgnlTTl one can get the Keldysh component Green function: 




(74) 



which is nothing but the FDT Eanl32lin any equihbrium system. Obviously, it is quite different than iG^ in Eanl32l 
for the non-equiUbrium system. 

Very similar expressions for the greater, the lesser, time ordered anomalous Green functions and the Keldysh 
component can be written down by changing the subscript n from Eqn. [7T] to Eqnl74] to the subscript to a. Two 
exciton correlation functions can be computed using the Wick theorems. Their analogies and differences with those 
of the excitons in non-equilibrium systems calculated in Sec. VIII can be addressed similarly. 

Because the retarded and advanced Green functions are the same as those in the non-equilibrium case, so the 
normal and anomalous spectral weights are the same as those in the non-equilibrium case listed in Eqnl52| 1541 and 
drawn in the Figs. 5 and 7. However, the distribution functions are quite different. For the equilibrium case in the 
Fig. 2, the distribution functions are just multiply Figs. 5 and 7 by sgn{uj). While for the non-equilibrium case in the 
Fig.l, the distribution functions are shown in the Fig. 6, 8 and 9. 



In this paper, by using the Heisenberg-Langevin equations in the context of the input-output formalism, we calcu- 
lated various one and two exciton correlation functions and also explored their relations to various photon correlation 
functions. We evaluated both the normal and abnormal spectral weights which lead to the excitation spectra of 
the non-equilibrium exciton superfluids. We also studied both the normal and abnormal Kelydsh component Green 
functions which lead to the distribution functions of the non-equilibrium exciton superfluids. These relations brought 
out the important connections between the properties of exciton superfluid systems to be probed by transport ea- 
surements and the various quantum optical photoluminescence measurements on emitted photons such as the angle 
resolved photon power spectrum, phase sensitive homodyne experiments and HanburyBrown-Twiss type of experi- 
ments. 

By the imaginary time formalism, we calculated the exciton correlations functions of an equilibrium exciton super- 
fluid subject to a photon dissipation bath. It was demonstrated explicitly that Gn/a^^n/a ^"^^ spectral weights 
Pn/a are the same for the two systems, but the G^^^, G^^^, G^^^ and the distribution functions, namely, the Keldysh 
component Green functions G^y^ are different. The analogies and differences between the two kinds of systems can 
be summarized in the Fig.l and Fig. 2. In the equilibrium system in Fig. 2, one measures Q"^ experimentally which is 
directly related to the spectral weight p by the FDT EqnlTT] ( see also^ ). The conventional theoretical tool is the 
imaginary time path integral method which can be analytically continued to the real time as shown in the Sec. IX. 
While in an non-equilibrium system in Fig.l, in the present context of non-equilibrium exciton superfluids, one also 
measures G^^ for emitting photons which is related to the corresponding G^ for the excitons by the input-output 
relation EqulSO] derived in Section VI. Unfortunately, due to the lack of the FDT, the G^ for the excitons is not 
simply related to the spectral weight p, its relation is non-trivial and need to be studied in details. For the given 
initial conditions EqnlSl this relation has been worked out by the Heisenberg equation of motions and was discussed 
explicitly in Sec. IV- VII. This important feature should hold in any non-equilibrium quantum open system. 

The input-output formalism can only be used to study the dynamic scattering process in Fig.l starting from the 
input initial state Eqn|5] While the imaginary time path integral method in the section IX can only be used to study 
the equilibrium case in Fig. 2 based on the exact ground state Ean l63l It was known that the correlation functions in 
both equilibrium and non-equilibrium quantum open systems could also be calculated by the Keldysh Green function 
method. The Kelydsh formalism in either Canonical Quantization language or path-integral language can be used to 
study both the dynamic scattering process in Fig.l starting from the input initial state Eqnl5]and the the equilibrium 
case in Fig. 2 based on the exact ground state Eqnl63l In a non-equilibrium case, because the G^^^ and G^^^ are 
independent of each other, one need both the forward and the backward paths in the real-time Keldysh contour. 
While in an equilibrium case, they are related by Fluctuation and dissipation theorem (FDT), so the forward and 
backward path can be reduced to just the forward path. In the Kelydsh formalism, one can find G^^^, G^^^, therefore 
the spectral weight — i{G^^^ — G^y^) — i(G^y^ — Gj^y^) which leads to excitation spectrum of the interacting 
system. Most importantly, one can also determine the Keldysh component G^y^ = ^n/a ^n/a ''^hich leads to the 
distribution function. It is the distribution function which lead to the difference between the equilibrium and the non- 
equilibrium system. Then one can determine G^, and G^, respectively and compare with those achieved from the 
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input-output formalism for non-equilibrium system and the imaginary time path integral method for the equilibrium 
system. Two-exciton correlation functions can also be evaluated by the Keldysh formalism by the counter-ordered 
Wick theorem. As shown in the Sect. IV- VIII, the procedures using the input-output formalism is just opposite to those 
using the Kelydsh formalism: one determine G^^^ and G^^^ first, then the spectral weight p„/a — i{G^^^^ — ^n/a)' 
then determine G^^^, G^^^, G^^^ and G^j^. The connections between the Keldysh Green function formalism and the 
input-output formalism will be explored in details in a future publication^^. 

The various one and two exciton correlation functions can be used to calculate the transport properties of the 
excitons such as the Coulomb drag and counterflow experiments. They are also needed to calculate the superfluid 
density and the critical velocity of the dissipative superfluids. Indeed, as pointed out in the introduction, the transport 
properties can provide more direct evidences on superfluids than the spectroscopic properties. The input-output 
relations can be used to establish the connections between the transport measurements on excitons and the PL 
measurements. In a future publication, using all the one and two exciton correlation functions derived in this paper, 
we will compute the Coulomb drags, the counterflow resistances, superfluid densities and critical velocities on both 
equilibrium dissipative superfluids in the Fig. 2 and the non-equilibrium open system in the Fig.l. 

As pointed out in the introduction, in this paper, we ignored the spins of the electrons and holes, therefore also 
the polarization of emitted photons. In fact, the electrons in the conduction band carry spin s — ±1/2. Due to the 
spin-orbit coupling, the holes in the valence band carry total spin J — 3/2. The quantum well confinement potential 
in the EHBL splits the energy of the light hole with m; — ±1/2 above that of the heavy hole = ±3/321. So we 
will neglect the higher energy light-hole, only consider the lower energy heavy- hole with rrih — ±3/2. One electron 
with s = ±1/2 and one heavy-hole with nih — ±3/2 can bind into four kinds of heavy-hole excitons (±1/2, ±3/2). 
When coupled to photons, the selection rule of the photon angular momentum is J = s + m^, so the 4 kinds of 
heavy-hole excitons can be grouped into the bight excitons J = ±1 which couple to the one photon process with the 
polarization cr = ± and the dark excitons J = ±2 which do not couple to the one photon process. However, the spins 
of electrons and holes inside an exciton will relax ( or de-cohere ) due to the spin-orbit couplings in the conduction 
band and valence band respectively. These individual electron and hole spin flips lead to the conversion from the 
dark excitons to the bright excitons. The simultaneous flipping of both electrons and holes is due to the Coulomb 
exchange interaction, it leads to a spin flip between J = 1 and J = — 1 inside the bright excitons. In order to make 
quantitative connection between theoretical predictions with the experiments, it would be important to study how 
the results achieved ini^iii and in this paper for spinless excitons will be modified after taking into account of both 
bright and dark excitons, also the effects of a trap, disorder and finite temperatures. 
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